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Using non-local chiral quark model with simple pole ansatz for mass dependence on momentum 
and non-local currents satisfying Ward-Takahashi identities, we calculate pion to photon transition 
distribution amplitudes and relevant form factors. For vector amplitude we recover correct normal- 
ization fixed by the axial anomaly. We find that due to the non-locality in the vector current, the 
value of axial form factor at zero momentum transfer is lowered with respect to the vector form fac- 
tor. Such behaviour - consistent with experiment - is not seen in the local models. Where possible 
we compare our results to the experimental data. 



I. INTRODUCTION 

Ov 

Transition Distribution Amplitudes (TDA) were firstly introduced in ref. [l| as objects parametrizing soft part 
of the amplitudes for hadron-antihadron annihilation HH — > 7*7 or backward Compton scattering j*H — > Hj, 
where the hard scale is provided by high virtuality of one of the photons. In some sense TDA's are hybrids of the 
ordinary Distribution Amplitudes (DA) and Generalized Parton Distributions (GPD) (see [2| for a review). However, 
if one restricts oneself to the mesonic case, they are more similar to GPD's - the difference is that we deal with 
— ^ , matrix elements which are non-diagonal not only in momenta but in the physical states as well. Nevertheless, from 
kinematical point of view, both are almost identical and therefore we can use similar variables as skewedness for 
example. In practice we consider two kinds of TDA's: vector and axial, depending on nature of the bilocal quark- 
antiquark operator sandwiched between photon and meson states. For experimental issues of Transition Distribution 
Amplitudes see ref. [3j]. 

In order to avoid complexity of bound state physics we limit ourselves to pions only. Pions are Goldstone bosons 
Q_i' of broken SU(2) chiral symmetry and their properties are to large extent determined by the symmetry (breaking) 
alone rather than by the complex phenomenon of confinement. After such simplification, we can use the non-local 
semi-bosonized Nambu-Jona-Lasinio (NJL) model to get an insight into TDA's. On the other hand, since Transition 
Distribution Amplitudes are related to anomalous diagrams, they can serve as a demanding tester of the non-local 
CN 1 models. We shall come back to this point later in this paper. 

First estimates of TDA's were made in refs. The first one IJ] was based on general QCD symmetries 

supported by simple quark model, while the second [g] and the third one [3p made use of Spectral Quark Model (SQM) 
and Pauli-Villars regulated NJL, respectively. The calculations in SQM respect all QCD symmetries by definition 
(i.e. Lorentz invariance, Ward identities, anomalies etc.). In Pauli-Villars regulated NJL model the presence of the 
finite regulator Q gives normalization for the vector amplitude which is not consistent with axial anomaly. In Q we 
calculated TDA's numerically in the non-local chiral quark model, however, since we have used "naive" currents we 
also broke the normalization condition of vector TDA. 
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In the present paper we extend our calculations from ref. [8| in twofold way: we use the full non-local vertices 
and obtain not only numerical but also analytical results for some of the TDA's. As a result we recover correct 
normalization conditions fixed by axial anomaly and give predictions for axial form factor which is not restricted by 
anomaly. Our results are here in qualitative agreement with experimental data. Since we revised our old calculations 
using different method we traced the mistake in the vector amplitude, which however does not change qualitative 
?3 , results. 

The paper is organized as follows. In Section[ll]we review the non-local chiral quark model and give the prescription 
for modified currents. Next, in Section IIII1 we recall the definitions of TDA's and relevant sum rules. Section [IV] 
contains our results with special emphasis on the role of the new pieces coming from the non-local parts of the currents. 
In Section [V] we investigate the relevant form factors and finally we summarize our results in Section FVll Technical 
details are given in the Appendices. 
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II. NON-LOCAL CO VARIANT CHIRAL QUARK MODEL 



In the following we use the non-local semi-bosonized Namu-Jona-Lasinio model. It is based on the quark-pion 
interaction in the following form [§| 



d A kdH 



Smt = J -^^-i>(k)^M{k)U^(k- l)^MiT)m, (1) 

where M (k) is dynamical quark mass appearing due to the spontaneous chiral symmetry breaking. Meson field J7 75 
is given in terms of the pion field as 

[/T5( ;E ) =exp |^_ r a 7r a (a;)75 J ) (2) 

where F ff = 93MeV is the pion decay constant. 
One defines mass dependence on momentum as 

M(k)=MF 2 (k) (3) 

where form factor F (k) should vanish for k — ► oo and is chosen to satisfy F (0) = 1. Expression for F (k) was 
obtained analytically in Euclidean space from the instanton model of the QCD vacuum and it is highly non trivial 
0| . Therefore here we use the Minkowski form proposed in [l(J 

/ -A 2 \" 

which reproduces reasonably well the instanton result when continued to Euclidean space. However, it can be also used 
directly in the Minkowski space. Integer parameter n defines a family of models and allows to analyze a dependence of 
our results on the shape of F (k). It was argued in ref. [llj that F (k) should vanish exponentially for large momenta 
in order to be consistent with OPE. This is another reason for introducing n parameter, which allows to control high 
momentum behavior. The method of fixing value of A„ is described below. 

It is a well known fact that due to the mass dependence on momentum the standard vector and axial currents are 
not conserved and Ward-Takahashi (WT) identities are not satisfied. Although it was argued in ref. [1] that this 
violation is not very large and almost does not affect shapes of TDA's, it certainly affects overall normalizations. It 
is therefore necessary to resolve the problem of normalization if one wants to use TDA's calculated in the non-local 
model for phenomenological estimates. 

There are several ways of constructing conserved currents in presence of non-local interactions see refs. [H. fl^. fl3. 

EE 

fl6| . However, it should be recalled at this point that none of them is unique, because current conservation fixes 
only the longitudinal part of a given vertex while the transverse one has to be modeled. In this paper we use the 
simplest "minimal" vertices satisfying WT identities. For vector current j M = {[ij^ip we replace 7 M by the following 
non-local vertex 

r»(k,p)=<f+gi>(k,p), (5) 

where the non-local addition reads 

9"(k >P ) = -^±^(M(k)-M(p)). (6) 

This form of non-locality does not introduce singularities as required by general prosperities of vector vertices 
Axial current j£ = ^7^75^ is made conserved by replacing 7^75 by 

r%(k,p) = j»j 5 +gg(k,p), (7) 

with 

& (k,p) = (M (k) + M (p)) 75. (8) 

(p-k) 



In contrary to the vector vertex, the axial one contains physical singularity corresponding to the pion. 



M = 225 MeV 


n = 1 


A = 1641 MeV 


n = 5 


A = 3823 MeV 


M = 350 MeV 


n = 1 


A = 836 MeV 


n = 5 


A = 1970 MeV 


M = 400 MeV 


n = 1 


A = 721 MeV 


n = 5 


A = 1704 MeV 



Table I: Numerical values of the model parameters obtained using Birse-Bowler formula (J9j) for the pion decay constant F^. 



For a model given by n we fix A„ using the Birse-Bowler |15l | formula for the pion decay constant F n : 



N c f°° jl2 M 2 (k E )-klM(k E )M'(k E ) + kl 1 M'(k E y 



4^ 2 J + M 2 (fc E )) 2 

with F n = 93 MeV. Using (@]) this expression can be calculated analytically - see the Appendix of jjjj ■ In Table [J we 
present A„ for given values of constituent quark mass M andpower n. Let us recall at this point that there exists 
another expressions for F„, namely the Pagels-Stokar formula [12] • As shown explicitly in e.g. [l8| it accommodates 
a non-local pion-quark interaction, however it corresponds to the naive axial current j§ which, as mentioned above, 
does not exhibit PC AC. On the contrary the Birse-Bowler formula takes into account the non-local axial current 
(more precisely its derivative, which is determined unambiguously). See ref. [l8| for discussion. 

It is important to note that A„ does not correspond to the QCD scale characteristic for the present model. Indeed 
function F(k) does not change much for different (n,A n ) except for the high momentum part. Therefore, as pointed 
out in [1, the precise definition of the scale can be done only within QCD, while within an effective model one 
can only estimate the order of magnitude. For the instanton model the characteristic scale is about 600 MeV [9|]. A 
somewhat more detailed discussion of this issue can be found in ref. (2(J • 

The non-local model with momentum dependent quark mass given by (J4j) was applied in the past to several low 
energy quantities (pion DA (lol |. two pion generalized DA [jj^], quark and gluon condensates [13]), using however 
naive vector and axial currents. Later in ref. fl^ |. pion DA with conserved axial current was estimated using PCAC. 
Recently a set of photon DA's up to twist-4 was obtained using currents satisfying WT identities [l7| . 

Similar models with however different form of F (k) in were also considered in the literature. In ref. [HI 
exponential form of F (k) was used. In refs. pll. [2!. 123. [2!. I25L l26| various distribution amplitudes and correlators 
have been studied in non-local models similar to ours. Pion and kaon DA were calculated also in 1271. 



III. KINEMATICS AND DEFINITIONS 



Relevant kinematics is very close to the one used in GPD formalism. We consider pion with momentum Pi and 
real photon carrying momentum P^. We shall work in the chiral limit, therefore both P 2 — and P 2 = 0. We define 
momentum transfer = P£ — -Pf and momentum transfer squared t — q 2 . Two light-like directions are defined by 
null vectors n — (1,0,0, —1) and h = (1, 0, 0, 1). Decomposition of any vector i> M in this basis reads 

«"=„+_+«-_+„£, (10) 

while the scalar product 

1 , _ 1 _ , _ _, 
u ■ v — —u v + —u v — Ut ■ Vt, (11) 

where arrows denote Euclidean two-vectors. Introducing average momentum p = h (Pi + P2) we can define so called 
skewedness variable 



(12) 



For massless pions — 1 < £ < 1. Switching to the frame where the average momentum does not have transverse part, 
i.e. 



+ iV p 2 n" 
P =P T VT 



we can write the following parametrizations of the pion and photon momenta 



(13) 



(14) 



r 1 



s j-n^ , n 2 1 „ 

(i-e)P + T + (i + o^ T + ^. 



Notice that p 2 — — i/4. We denote photon polarization vector by e. It satisfies condition e • P 2 = 0. 
We can now define TDA's. Vector Transition Distribution Amplitude (VTDA) V(X,£, t) is defined as 



dX 
2^ 



7(^2,e) 



,/ j _ V) r»u 



• </'i>)= 2 ^ + e^et Pa q V (A, £, t) . 



In the axial channel the Axial Transition Distribution Amplitude (ATDA) A(X, £, t) is defined as: 



(15) 



(16) 



i\x P + 

2tt 



7(^2,e) 



7T+ (A) 



if 



2 2V2F w sign(0^ fX+j 
2£ 



(17) 



Notice that we have explicitly written down the contribution containing massless pole 1/t coming from the fact that 
the pion can couple directly to the axial current. This part is connected with the pion DA (/>„ (u). Dots stand for the 
remaining parts which vanish when contracted with n M Q. 
There exist the following sum rules for TDA's 



J I A(X,U) "V \ F A {t) 



(18) 



where Fy (t) and Fa (t) are vector and axial form factors respectively. Quantity Fy A (t) — Fy.A (t) /m^ is introduced 
because we work in the chiral limit (m ff is the mass of the pion) . Moreover the transition form factor for the process 
7t° — > 7*7 is related to the vector form factor by the formula 



F W1 (t) = V2F V (t) . 



(19) 



The value of F ni for t — is fixed by the axial anomaly and equals 



(0) 



1 



0.272 GeV -1 . 



(20) 



There is no such a constraint for the axial form factor. See section [V] for the discussion. 



IV. TDA'S IN COVARIANT NON-LOCAL MODEL 



In this section we outline some points of the calculations within the model introduced in the previous Section and 
show our results. More details are given in the Appendix lAl and 151 



A. Vector TDA 



In the vector channel direct calculation of the matrix elements gives 



d { -~n) Vui-n 



2 \2 



7T+ (Pi) 



V 2< } ' Vc (Q d M^ (X, £, t) + Q U M V V £, t))et, (21) 



where N c is number of colors, Q u , Qd are charges of the pertinent quarks and 

Tr {5 (fc + P 2 ) r^ 1 (fc + P 2 , k + Pi) 5 (fc + Pi) T 5 (k + P u k)S (k) V (k + P x , fe)} 

with 

SM = j^¥r (22) 

In the formula above quark-pion coupling reads 

T 5 (k,p) = F(k)-f 5 F{p). (23) 

Notice that we have used the non-local vertex |(5]) also in the bilocal currents. This is necessary if one wants to 
maintain WT identities. 

Amplitude My" , due to the natural splitting of the vector vertices to local and non-local parts, can be divided into 
four pieces 

M v v = M v v - (0) + M v v ' (1) + My"' (2) + My' (3) . (24) 
To list them we introduce shorthand notation 

dk = -0^S(k + -(X-l)p+). (25) 



Then 



AC (0) = j dk r Ii{S{k + P 2 )^S{k + P l )T b (k + P 1 ,k)S(k) 1 u } (26) 
corresponds to the quantities calculated in @|. The new pieces are 

M^ v,(1) = [ dkTr{S(k + P 2 ) 1 > i S(k + P 1 )r 5 {k + P 1 ,k)S(k)g u (k + P 1 ,k)}, (27) 



My' (2) = J dk Tr {S (k + P 2 ) g"(k + P 2 ,k + P 1 )S(k + P x ) T 5 (k + P u k)S (k) 7"} . (28) 

The remaining part of My' is zero, simply due to vanishing of the Dirac trace, M 1 ^' ^ = 0. 

Above decomposition, after projecting on the proper tensor structures, leads straightforwardly to the expression 
for VTDA 

V (X, C, t) = V<® (X, C, i) + V<» (X, t) + (A, t) . (29) 

The V^+V^ part is the addition that is required in order to recover correct normalization. The explicit expressions, 
together with some details needed to calculate analytically are given in the Appendix [Aj 

Our results are shown in fig. [TJ We plot also separately the local part and the addition + itself. We 
compare the new result with the local version of our model calculated in [8]. Qualitative behavior with M, £, n and t 
is the same as in [8j therefore we do not discuss this further. However now the condition coming from axial anomaly 

V{X,tt) = ^ (30) 
is satisfied automatically for any value of n and A (which does not have to be necessarily equal A„). 
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Figure 1: Vector TDA for M = 350 MeV, t = -0.1 GeV 2 and n = 1. a) Solid line representing full VTDA for £ = 0.5 is the 
sum of the part coming from the local part of the vector vertex (dotted) and the non-local one (dashed); b) comparison of the 
results of this paper for £ — 0.5 (dashed), £ = —0.5 (dotted) and the local version of the model M (k) = M for £ = ±0.5. 



Axial TDA 



For the matrix element of the axial operator we have 



/ 



where 



d — — n u 



7T+ (Pi) 



y/2eMN c 



(Q d M^ (X, e, t) - Q U M% (-X, £, t)) el, (31) 



(27T) 4 



<5(fc+-(A:-l)39 + ) x 



Tr {5 (*; + p 2 ) r£ (fe + p 2 , k + 1\) s(k + Pi) r 5 {k + P!,k)s (k) r» (k + p u k)} . 



This expression looks very similar to (|2T]l , however here full axial vertex Tg given in ([7]) appears instead of the vector 



one. 



Exploring this expression we find that it has again four parts, similarly to the vector case. However only two of 
them contribute to the ATDA, namely: 



and 



M 



M 



f", (0) 



= / dkTi{s(k + p 2 )^ l5 s(k + p 1 )r 5 (k + p 1 ,k)s(k)Y} 



= / dkT r {S(k + P 2 ) 1 » l5 S(k + P 1 )r 5 (k + P 1 ,k)S(k)g»(k + P 1 ,k)} 



(32) 



(33) 



The remaining terms are connected with pion DA (or they are gauge artifacts). This can be seen by noting that 

gg (k + P 2 , k + Pi) = ^ (M (k + Pi) + M (k + Pi)) 75 (34) 

and comparing with definition l|17p . Therefore we have 

A (X, t) = AW (X, £, t) + (X, t) , (35) 

where again A^ is the old part already calculated in [£]. All further details concerning calculations are relegated to 
Appendix [Bj 

The results are shown in fig. [2j It turns out that the integral over dX of A^ (X, £,t) part is negative, what shifts 
the value of the axial form factor towards the experimental value. For negative values of £ the addition A^ (X, £, t) 




Figure 2: Axial TDA for M = 350 MeV, t = -0.1 GeV 2 and n = 1. a) Solid line representing full ATDA for £ = 0.5 is the sum 
of the part coming from the local part of the vector vertex (dotted) and the non-local one (dashed). Dash-dotted curve is the 
results for the local model M (k) = M. b) The same for £ = —0.5; c) results for positive £ = 0.5, n = 1, t = —0.1 GeV 2 and 
different values of M; d) the same for negative £ = —0.5. 



has in a sense unexpected shape, what quite drastically changes the shape of full ATDA for £ < 0. This effect is 
weaker for lower constituent quark masses as the addition is proportional to the third power of M/A n . We have 
checked that for any £ (both positive and negative) 

J dX A™ (X, |e| , t) = j dX (X, - |£| , t) = / (t) < 0, (36) 

where / is some function of t only. Equation (f36|) is a special case of the polynomiality condition, which is satisfied 
in our model, both for VTDA and ATDA. 



V. FORM FACTORS 



As already mentioned in Section llV Al vector and transition form factors at zero momentum transfer are reproduced 
correctly, as given by axial anomaly. Full result for transition form factor multiplied by momentum transfer is shown 
in fig. 03 We make a comparison with CELLO [H], CLEO and some new BaBar data [3(3] for several values of 
model parameters. Since our approach is reliable at low momenta transfer we limited ourselves to the region of a few 
GeV 2 . We do not use the standard dipol parametrization for the transition form factor [29J as a reference, because 
recent BaBar data seem to reveal different behavior (see also [3l|, HH). We find that model parameters M = 300 MeV 
and n = 1 fits the data quite well, while low constituent quark masses lie much below. However, we underline that 
momentum transfer around a few GeV 2 can be too large to be treated within the considered approach. 



In the case of the axial form factor, we find that due to the contribution of A^- 1 ' its value at zero momentum transfer 
Fa (0) is lower than for the vector one. The values of Fa (0) and the relevant ratios Fa (0) jFy (0) for several model 
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Figure 3: Pion to photon transition form factor times momentum transfer for several model parameters versus CELLO, CLEO 
and BaBar data in the — 8 GeV 2 range. 



M [MeV] 


n 


F A (0) 


F A (0) /F v (0) 


225 


1 


0.0217 


0.80 


350 


1 


0.0168 


0.62 


350 


5 


0.0163 


0.60 


400 


1 


0.0161 


0.60 


400 


5 


0.0152 


0.56 



Table II: Numerical values for axial form factor at zero momentum transfer Fa (0) for different model parameters. 



parameters are presented in Table HH We notice that when n is increasing, the value of Fa (0) is decreasing very 
slowly. Therefore we conclude that Fa (0) is quite robust as far as model parameters are concerned. Experimental 
value given by PDG is 

F^ p (0) = 0.0115 ± 0.0005 (37) 

and 

(F A (0)/F v (0)) cxp = 0.7t° 6 2 . (38) 

Although our predictions are much better then the ones obtained in local models (with Fa (0) = Fy (0)) 0, [H, 0, @|, it 
still does not agree with experimental value. Notice however that also prediction for Fy (0) provided by conservation 
of vector current and axial anomaly 

F v (0) sa 0.027 (39) 

overshoots the experimental value (PDG) 

Fy p (0) = 0.017 ±0.008. (40) 

We recall also that our calculations were done in the chiral limit. The effects due to the finite pion mass were discussed 
in OH within local NJL model and turn out to be small. 



VI. SUMMARY 



In the present paper w have used the instanton motivated non-local chiral quark model to compute pion-to-photon 
transition distributions amplitudes in vector and axial channel defined in eqs. (fl6|l and (fl7|) . Our approach follows 
closely the approach of ref. [1] with one important modification. In ref.pj we have used the momentum dependent 
constituent quark mass and the nonlocal pion-quark interaction, however, both vector and axial currents were not 
modified to satisfy Ward-Takahashi identities. We have argued in @ that such an approach - although inconsistent 



- gives reasonable shape for the TDA's missing, however, normalization. Indeed, our present findings confirm this 
conclusion. 

In order to satisfy WT identities the "naive" QCD currents require modifications if the quark mass depends on 
momentum. Such modifications are not unique because WT identities do not fix the longitudinal part of the current. 
In the present paper we have used the simplest possible modifications of the vector and axial currents given by eqs.((5]) 
and 10. 

More importantly, we have used the nonlocal vertex not only for the photon coupling to the quark loop but also for 
the vector operator defining the TDA. This modification is not obvious for the following reason. Factorization of the 
physical process into hard and soft parts is done in QCD in terms of the operator product expansion and the relevant 
operators entering the definition of the soft part are therefore indeed the "naive" QCD operators. In our language 
this would correspond to neglecting M^f' ™ of eq. f28|) . We have checked that this would lead to the violation of the 
normalization condition for VTDA. In the case of ATDA the non-local part of the axial current does not contribute, 
since it enters the piece proportional to the pion DA - see eq.JTTJ) - but not the ATDA itself. 

Our findings can be summarized as follows. The normalization of TDA's in the vector channel - which is fixed 
by the axial anomaly - turns out to be correctly reproduced in our approach. At the same time normalization of 
ATDA's is lowered due to the non-local term in the vector current. As a result the two normalizations are no longer 
equal, as observed experimentally. The shapes of the TDA's are similar to the ones obtained in the local model and 
in the non-local model with naive currents. In the case of VTDA's we observe, however, small difference between 
distributions with positive and negative £ what was previously not the case. ATDA's acquire negative contribution 
from the non-local piece which makes the differences between positive and negative £'s even more profound than in 
the previously considered models. We have also investigated dependence of TDA's on the choice of constituent quark 
mass M and the shape of the cutoff function F(k). Dependence on the power n entering F(k) is rather weak, whereas 
the dependence on M is effectively important only for ATDA for negative £. 

We have also calculated the corresponding form factors. Axial form factor is quite robust with respect to the shape 
of the mass dependence on momentum. Pion-to-photon transition form factor, which is directly related to the vector 
form factor, can be compared with experimental data. We compare it with CELLO, CLEO and recent BaBar data 
in low energy regime in fig. (J3j) and conclude that our results reproduce quite well its t dependence. 

Let us finish by a remark that we have worked in the chiral limit where — 0. The effects of adding small current 
masses are expected to be small, however, only further study can show whether they will be able to correct small 
deviations from experimental values of Fv,a{0). Finally the modifications of the nonlocalities that are still possible 
as far as longitudinal part of the currents is concerned deserve further investigation. 
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Appendix A: VECTOR TDA 

In the beginning we introduce a shorthand notation in order to make the formulae more compact: following 
abbreviation 

M (k) = M k (Al) 
and similarly for other fc-dependent quantities. We recall that Mk = M Fj*. For the inverse scalar propagators we use 

D k = k 2 - Ml (A2) 
Let us now turn to VTDA. We obtained the following expressions for subsequent pieces giving contribution to 

v(x,s,t) 

V® (X, £, t) = 16t MN cP + (Q d I {l) (X, e, t) + QuT {l) (-X, & i)) (A3) 

where 

I<°> (X, £, t) = f dk FkF ^ x J Fl ( Mfc+Pl - M k+P2 ) + F 2 (M k - M k+P2 ) + M k } , ( A4) 



In the above formulae 



J -Dfc+Pi -D/c+p 2 D k (k- P 2 ) 

IV (A, t t) = -l dk (-^ 2 °T)F k F k+pA M, +P -M k+P ^ 
J D k+Pl D k+ p 2 D k (k- q) 



and 0t is an angle between kr and qr in the transverse plane. Integration measure dk was introduced in eq. (|25|) . 

As an example we briefly explain how to calculate 1^. Following refs. [lCJ] we introduce dimensionless variables 
n = k/A, Pi = Pi/A, P 2 = P 2 /A, r = M/A and use 

"i,2 = (« + Pia) 2 - 1 + (A9) 
u 3 = K 2 -l + ze. (A10) 

Then using (J4j) we have 

^ (A, 6* = -T3 / d« T . - ^ 1 2 3 1 2 j , All 

where 

G(z) = z 4n+1 + z 4 "-r 2 . (A12) 

We can obviously write 

4n+l 

G(z) = J] (2-77,), (A13) 

8=1 

where rfr-S are the roots of G (z) — to be obtained numerically. Their properties were discussed in [JjJ] and [l7| . 
Then using decomposition into proper fractions we obtain 

4n+l 

Z« (A, = £ fifjf^ijklSl (A, & t) , (A14) 

where 

4n+l 1 

; -=nra ,AI5) 

t>«s = lf"lf <fi ('ff - >)?) (A16) 

and 

m^=f & pi 1r 2 ^u — 7- (Ai7) 

Next we assume that this expression can be continued to the Euclidean space - this is equivalent in Minkowski space 
to deformation of the dn~ integration contour as described in [Tj], [l?]]. This choice of the contour assures that the 
results are real and analytical in A„. The remaining parts of integrals can be transformed to a sum of 



the integrals similarly as in the case It is convenient to start with the dn~ integration using residue theorem. 

Positions of the poles in the k~ complex plane depend on the kinematical region and on the sign of £, however the 
pole connected with k ■ P2 does not give contribution as it should be. Physical support — 1 < X < 1 splits into three 
pieces: — 1 < X < — |£|, — |£| < X < |£| and |£| < X < 1. Integrals (X,£,t) are nonzero in the last two, while 
XW (—X,£,t) in the first two intervals. 

In order to write down the results it is useful to introduce the following notation 

g\ A) = (4 + 1) (i - - fa ■ q T (x + 1) - (x - 1) (1 - 3£ - x (1 + 0)P 2 + (i - vu (ais) 
g\ B) = (4 + 1) (i - - «t • <!t {x-i)-{x- if (1 + of + (i - (Ai9) 

^ = (4 + 1) £ + * • <7T - (X 2 - 1) £p 2 + , ( A20) 

= (4 + 1) £ + (* - 1) «t • q T -{x- if CP 2 + fa, (A21) 

d (AB) = ^ ^2 + ^ _ 2X Ht . ^ + 2 (x 2 - 1) CP 2 + (X-S)n-(X + (A22) 
djf C) = - (1 + (4 + 1) - (1 - X) K T ■ q T + (1 - (X - iff + (X + % - (X - 1) % , (A23) 

4? B) = (! - (4 + 1) + (i - x) k t ■ q T + (i + o (x - iff + (x - o m - (x - 1) (A24) 

Ajfe = -a/y» (A25) 
where p — p/ A is dimensionless average momentum. Then we have 
• |£| < X < 1 

z»> K| , ,) . r»> (X, - 151.') = ( x ~ l ) E /.AA / £S ^y/^ <A26) 

i,j,fe \ I a ki kj 

4n+l „ ,2 2-2/j 

X« (X, Id , t) = I« (X, - |f| , t ) " (X - l) 2 2 ^/^A / ^ (^S^j (A27) 

^ i,j,k J y Z7T ) 9k a ki a kj 



1 ■ 4n+l . , 2 2-2/j 

K») (X, |C| , t ) = *»> (X, - |f| ,t) = 1^ (X - if £ ^/^A / ^ ^ggg^) (A28) 



4n+l „ ,9 2-2 

v r f f f / t -4 sin 1 

ki<*<ki 

4,1+1 /• j2 ., p„ i n , .. ET i „3n„4n„ra 



H (X, \$,t) (X + 0^ fifih y r— 3 jab) (_AAC)\ (A29) 

i,j,k y > u ij y u ki j 

1 - Id , *) = - g a^T ( X - ^) A 7 d (AB) d (CB) ( A3 °) 



' d Kt —KTp sin &t 



i,j 7 k v / a A;j 



Wlf|,t) = -i^ (* + a EA»/</iA /^^tS^ (A33) 



(2») 3 (-4f c| ) 



i w W-KI.0-5^:('-O , E*.AAA f < A34 > 

Integral Zw can be calculated completely analytically, however we do not present the result here to make the 
paper more compact. It should be pointed out that the analytic formulae contain large sums of complicated complex 
expressions - therefore it is usually more efficient to perform all integrals numerically. 

Appendix B: AXIAL TDA 

In the axial channel we get 

A« {X,£,t) = -WiMN c p+ (Q d J® (X,£,t) - Q u J {i) (-X,£,t)) , (Bl) 

where 

JW (X, £,t)= [ dk ffcFfc+Pl 

x jjl [M fe+Pl - M k+P2 - 2Af fe ] - F 2 {M k+P2 - M fe ) + M fc + 2F 3 [M k+Pl - M k ] J , (B2) 



(x,z,t)= I dk FkFk+p > {Mk ~ Mk+P > } 



• P2D k+Pl D k+P . 2 D k 

x jVi [M fe+Pl Mfe - Mfc+P 2 Mfe - 2/c 2 - 2fc ■ p] + F 3 [M k+Pl M k+ p 2 - M k M k+P . 2 + M k M k+Pl 

-fc 2 +P 1 -P 2 ]|. (B3) 

Functions Fx, F2 are the same as in vector case while 



8£ (k T ■ qr) +2tk T - q T (X - 1) (2e +£-l)+t(e-l)(t(X- l) 2 (1 + - 4£fc 2 
F 3 = — i '- « i '-. (B4) 

2(1-0 i + O * 2 



Calculations proceed quite similarly as before, although they are more involved. However, in the case of 
additional difficulty arises if one wants to obtain result for any n. Let us discuss this further. 

Problematic part comes from the subterm of containing more then 4 powers of F(k). Using (Sal, {EEJ this 
part reads 



m , r 3 f - u^ul n (ul n - u?™) (F x + F 3 ) 

(X,U) = / dk 1 2 y 1 \\ ■ ) 3 j . B5 
A z J k ■ P 2 G (ux) G {u 2 ) G (u 3 ) u§ 



Notice that there is an additional u 3 in denominator which gives rise to multiple poles if n > 2. However if n = 1 

(«,.) =£±± nnh f dk w-y)(*+*> (B6) 

A ijk ■> K ■ P 2 (Ul - Vi) {U2 - Vj) (U3 - Vk) 

where 

/<= n 7^70 (B7) 

and 770 = 0, while 771 , ... , 775 remain solutions of G = 0. In order to avoid calculating residues of the multiple poles 
if n > 2, we choose to close the contour in such a way that multiple pole does not lie inside. The price we pay is 
that relevant expressions become more complicated. Moreover, although the pole coming from k ■ Pi does not give 
contribution to it does contribute to and J^ 1 separately where = + J^p ■ This fact has to be 

taken into account when we decide to choose different contours for and . 

Now we are in a position to write the results. We use the same abbreviations as in the vector case, additionally we 
define: 

c ijk = —Fi {f3 ijk + lijtnfrg) - 2F 3 p ijk + F 2 a ljk + r/f^f rft (B8) 



X B ' C) = Fir 2 p ljk - F&jktV'W + F 3 r 2 (a ijk + (3 ijk ) + F 3 f3 ljk s^ B > c \ (B9) 



where 



t(A) = - (1 + 2Q 4 + (1 + rn kt ■ gr) (2X 1) + (X 1) (X (2j - 1) + l)p 2 (mQ) 



\p 2 



- (1 + 4 + (1 + Vi - k t ■ q T ) (X-1)-(X- 1) ((X - l) 2 (1 - + 2 (X + 0) ; 

s = YTl (B ! 1 ] 



t(B ) = (2i-l)K 2 T + (l + r, 3 +n T -q T ) (2X 1) - (X 1) ((2X - 1) (1 + Q + X fl f 

X — £ 

(-i + 04 + (i + ^ + «T-gT)(^-i)-(x-i)((x-i) 2 (i + e) + 2(x-e))p 2 

S < B > = — ^ ^— (B13) 

t <*> = 4 + (1 - 1} + (X - (B14) 

— 1 

s {C) = l + 7]k-2p 2 . (B15) 
We recall at this point that 770 = 0, therefore although sums below run from unity, notation like "dio" makes sense. 

• |£| <X< 1 

4n+l . , 2 

<7 (0) {X, 1(1 , t) = (X, - 1(1 , t) = J£ (X - 1) £ /i/i/* / ^ ^cwcb) ( B16 ) 

^ i,j,fc 7 \ A7T > a ki d kj 



{Fx + F 3 



. 3 4n+l „ ,2 

jP \z\ , t) = j« ( x, - lei , t) = £ ./•,./•,/, - ^ ) J |^ 



5 (-4? c) ) (-4f >)" sfMfMf ) HD" U w ) ^£ B) H fl,r 



(B17) 

jjp (x, if i , t) = j$> (x, -\z\,t) = £;(x- if Y: fififk / 7^3 mOTB) (B18) 

1 i,j,k ■' 9k d ki d kj 



\Z\<x<\Z\ 



(x,\^\,t) = -^(x + of^fif j f k I 



4n+l „ , 2 



(B19) 



4n+l . , 2 

i,j,k V ' u ij 



4n+1 "dV (^ 1+J F 3 )(^ + O rl+2 



7 (1) fX If I - -— V f f ft (n n n 4n - n 3n n 2n ) / ^ 1+f3 ^ At ^ (B21) 

^ ( ' lel ' } " ^ + h 3 { l 3 1 3 ] I ^ g mt B) (-4f (-4f 0" 



4 1} - ici ,*) = E ^ (« n - / (fl+F3)(x '°" +2 



t* Khh h " J J (^)%f) (_,#->) (-4D 



(B22) 



4n+l , 2 U) 

rf> (x, k, , t ) _ ( x + o" E WyA / ^ ,7 , (B23) 

v ; tfi a ij \ u ki J 



4,1+1 r ri 2 u> (B) 

4" (x, - m ,o - ( x + o- g ma / ^ - ^ KCB)) (B2 4, 

Notice, that some of the integrals are superficially divergent. Convergence is assured by the identity 
see [13] for the proof and discussion. 
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